Density-wave phases of dipolar fermions in a bilayer 



F. M. Marchetti^ and M. M. Parish^'^ 

' Departamento de Fisica Teorica de la Materia Condensada, 
Universidad Autonoma de Madrid, Madrid 28049, Spain 
^Cavendish Laboratory, J J Thomson Avenue, Cambridge, CBS OHE, United Kingdom 
'^London Centre for Nanotechnology, Gordon Street, London, WCIH OAH, United Kingdom 

(Dated: July 18, 2012) 

We investigate the phase diagram of dipolar fermions with aligned dipole moments in a two- 
dimensional (2D) bilayer. Using a version of the Singwi-Tosi-Land-Sjolander scheme recently 
adapted to dipolar fermions in a single layer [M. M. Parish and F. M. Marchetti, Phys. Rev. 
Lett. 108, 145304 (2012)], we determine the density-wave instabilities of the bilayer system within 
linear response theory. We find that the bilayer geometry can stabilize the collapse of the 2D dipolar 
Fermi gas with intralayer attraction to form a new density wave phase that has an orientation per- 
pendicular to the density wave expected for strong intralayer repulsion. We thus obtain a quantum 
phase transition between stripe phases that is driven by the interplay between strong correlations 
and the architecture of the low dimensional system. 



Density-wave phases such as stripes are apparently 
ubiquitous in nature. They are typically found in quasi- 
two-dimensional or layered materials where they 
manifest as periodic modulations of the electron density 
within the two-dimensional (2D) layers. Moreover, such 
stripes have been linked with high temperature supercon- 
ductivity jUini- However, despite their ubiquity and po- 
tential importance, their origins and behavior are still un- 
der debate. Indeed, a central question is whether stripes 
are driven by electron-electron repulsion or simply by the 
architecture of the underlying crystal structure [6 . 

One route to gaining insight into the problem is to 
study cleaner, more tunable analogues of these electron 
systems. Quantum degenerate Fermi gases with long- 
range dipolar interactions [H |8] provide just such a sys- 
tem in which to investigate density-wave phases. Such 
dipolar Fermi gases have recently been realized exper- 
imentally with both magnetic atoms [5] and polar di- 
atomic molecules [TUHIS- particular, ultracold polar 
molecules of ""^K ^''Rb have been confined to 2D lay- 
ers using an optical lattice [13], thus paving the way for 
exploring long-range interactions in low dimensional sys- 
tems. 

For a 2D gas of polar molecules, the dipole-dipole inter- 
actions can be controlled by aligning the dipole moments 
with an external electric field. For small dipole tilt angles 
6 with respect to the plane normal, the dipolar interac- 
tions are purely repulsive, while for 6 > 7r/4, the interac- 
tions acquire a significant attractive component such that 
the dipolar Fermi system is unstable towards collapse for 
sufficiently strong interactions [l4l - fT7] . Away from col- 
lapse, in the repulsive regime, previous theoretical work 
has predicted the existence of a stripe phase , even 

for the case where the dipolar interactions are isotropic 
{9 = 0) and the system must spontaneously break rota- 
tional symmetry [16]. Here we investigate the effect of 
the low dimensional architecture on density instabilities 
by considering dipolar fermions in a 2D bilayer geometry. 



We determine the phase diagram of the bilayer sys- 
tem within linear response theory, using a version of the 
Singwi-Tosi-Land-Sjolander (STLS) scheme [IH] recently 
developed in Ref. [16]. Based on this analysis, we show 
that the bilayer geometry can actually stabilize the col- 
lapse of the 2D Fermi gas to form a new density wave 
(Fig. [T]). However, in contrast to the stripes in the re- 
pulsive regime, this new stripe phase has density modu- 
lations along the direction of the dipole tilt (Fig. [2]) and 
can also be well described by a simplified STLS theory 
that involves exchange correlations only. Our work thus 
reveals a new quantum phase transition between two dif- 
ferent stripe modulations, where one phase is driven by 
strong repulsive correlations and the other is driven by 
the bilayer architecture. 

In the following, we consider the bilayer geometry 
shown in the insets of Fig. [2j Here, the dipole moments 
(of strength D) are aligned by an external electric field E 
lying in the x-z plane and at angle 6 with respect to the z 
direction. We parameterize the x-y in-plane momentum 
by polar coordinates q — {q,(t>), with (j) — correspond- 
ing to the direction x of the dipole tilt. The remaining 
system parameters are the bilayer distance d, and the 
Fermi wave vector kp — \/47rn {n is the density in each 
layer). For dipoles confined in a layer of width W, in the 
limit qW ^ 1, the effective 2D intralayer interaction can 
be written as |20) : 

vn{^)^Vo~2^TD^qaO,(^>) , (1) 

where £,{9, 0) = cos^ 9 — sin^ 9cos^ (f>, and Vq is the W- 
dependent short-ranged contact interaction. The confine- 
ment width W provides a natural cut-off for the quasi-2D 
system: A 1/W :$> kp. 

Likewise, in the limit W <ti d, we can write the inter- 
layer interaction as |5T]: 

vi2{q) ^ -2TTD^qe~'''^[^{9,(l3) +isin29cos(j)] . (2) 

Note that ior 9 ^ 0, this interaction is complex and sat- 
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isfies f2i(q) = i'i2(q) = Wi2(— q)- This arises from the 
fact that the interlayer interaction in real space is not 
invariant under the transformation r iH- — r. 

Assuming identical layers, one can parameterize the 
bilayer system using only three dimensionless quantities: 
The tilt angle 9, the bilayer distance kpd, and the in- 
teraction strength U = mD^kp /fi^ , with m being the 
fermion mass. The cut-off A and the contact interaction 
Vq should not be relevant since these do not affect the 
low energy behavior of dipolar fermions, and indeed the 
procedure we employ preserves this. 

We now turn to the linear response theory used to ana- 
lyze the inhomogeneous phases of the dipolar system. In 
the bilayer (and multilayers generally), the linear density 
response Srii to an external perturbing field Vf^* defines 
the density-density correlation function matrix Xij: 
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Sni{q,LLi) 
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where i, j are the layer indices. For a non-interacting gas, 
we clearly have Xij = (^ijll, where the non-interacting in- 
tralayer response function Tl{q, uj) can be evaluated ana- 
lytically [22] . Typically, one includes interactions via the 
Random Phase Approximation (RPA), where one uses a 
perturbing field that contains an effective potential due 
to the perturbed density: K,*^^* i-> Vj'^'^ + J2j VijSrij, with 
intralayer potential W22(q) = vii{q). However, as has 
been argued recently for the single layer case, RPA is 
never accurate for dipolar interactions, since it neglects 
exchange correlations |17l 118] which are important even 
in the long- wavelength limit [T5] . 

A straightforward and physically motivated way of in- 
corporating correlations beyond RPA is by means of local 
field factors Gij{q) (for an introduction to this method 
see, e.g., Ref. [13]). Here, the (inverse) response function 
now reads: 



X 



U{q,uj) 



«.,(q)[I-G.,(q)] . (4) 



Note that we clearly recover both RPA and the non- 
interacting case if we take, respectively, Gij = or 
Gij = 1. This response function can be related to 
the "layer-resolved" static structure factor Sij (q) by the 
fluctuation-dissipation theorem: 
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(5) 



In turn, we can approximate the local field factors using 
the STLS scheme IHj: 



Gy(q) 



dk q ■ k Vjj (k) 
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The response function Xij (and associated structure fac- 
tor Sij) can now be determined by solving Eqs. (|4|-([6| 
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FIG. 1: (Color online) Phase diagram for a dipolar Fermi gas 
in a bilayer at fixed interlayer distance, kpd = 2, as a function 
of 9 (see Fig. |2| and interaction U = mD^kF/tt^ ■ The liquid 
phase is superfluid (SF). The (green) open triangles [circles] 
set the boundary of the stripe phase oriented along = 
[(j) — n/2], derived from a self-consistent STLS calculation. 
The filled (green) square at dc — 0.75 and U ~ 15.65 is a 
quantum critical point beyond which there is a phase transi- 
tion between the two stripe phases. The (blue) open diamonds 
for the = stripe phase are instead determined including 
exchange correlations only (see text). These boundaries can 
be compared to the = 7r/2 stripe transition (dashed line) 
and the collapse instability (dashed-dotted line) for the single- 
layer case [THj . The shaded "bosonic" region is where the sys- 
tem can be described in terms of interlayer bosonic dimers. 
The (red) filled diamond and thick (red) line at 6 = 7r/2 in- 
dicate collapse in the bilayer. 



self-consistently. The STLS scheme has been heavily 
utilized for Coulomb interactions and it has proven to 
be very successful for describing the dielectric function 
of several strongly-correlated electron systems (see [23] 
and references therein). Following Ref. [16], we consider 
an improved version of the STLS scheme that has been 
adapted to the dipolar system. In essence, it ensures 
that our results are insensitive to A and Vq, by requiring 
that the intralayer correlations be dominated by Pauli 
exclusion at large wavelengths 2kp. 

For identical layers, we can assume that ^22 = S'n, 
^21 = SI2 (and similarly for the local field factors Gij). 
Note that the complex form of the interlayer potential pi) 
means that the interlayer factors 512 (q) and Gi2(q) are 
also complex. However, the symmetry Wi2(— q) — '^i2(q) 
is also preserved for both factors at each iteration step of 
our self-consistent scheme. This guarantees that physical 
quantities such as the "layer-resolved" pair correlation 
functions, gy(r) = 4j (■j/4(r)?/>Vo)-!/'j(O)'0i(i")), where 



1 + - 
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are always real, even when i ^ j. 

We determine the density instabilities of the bilayer 
system by analyzing the divergences of the static response 
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function matrix Xij {^.^ 0) • Specifically, we search for zeros 
of the largest inverse eigenvalue, 

= ^ - v,i[l - Gn] + |z;i2[l - Gi2]| . (8) 

A zero of x+'i^.i 0) at a critical wave vector signals an 
instability towards the formation of a density wave with 
period set by Qc. If the instability occurs for a specific 
direction ^, then the density-wave phase corresponds to 
a one-dimensional modulation (or stripe phase) of period 
27r/(7c oriented along <f>. In this way, we obtain the phase 
diagram plotted in Fig. [T] for kpd = 2. 

For tilt angles 9 < 6c — 0.75, we find a stripe phase 
along = 7r/2 that is of a similar nature to the one 
found in a single layer (dashed line of Fig.[T]). In particu- 
lar, it is driven by strong intralayer correlations induced 
by the repulsive part of Vn, as evidenced by the rela- 
tive insensitivity of Qc to the bilayer geometry and 9 (see 
Fig. [2]). However, the presence of the second layer can 
decrease the value of the critical interaction strength Uc 
for stripe formation, as one might expect from the form 
of Eq. ([S]). The attractive part of 1^12 (q) also ensures 
that the density waves along (p = 7r/2 in each layer are 
in phase. Similar results were found using the conserving 
Hartree-Fock (HF) approximation [18l|24], but for much 
smaller values of Uc, like in the single-layer case. The 
shift of Uc due to the other layer is relatively small for 
distance kpd — 2 (see Fig. [T] at small values of 9), but 
it can become substantial for smaller kpd since Eq. ([s]) 
depends exponentially on the bilayer distance. However, 
for smaller distances, we then encounter phases involving 
strong interlayer pairing p5H27] and the system would in- 
stead be better described in terms of interlayer bosonic 
dimers, as we discuss later. 

In the isotropic case {9 = 0), we find that the sys- 
tem spontaneously breaks rotational symmetry to form 
a stripe phase at [/ ~ 5.74, similarly to the single-layer 
case |16j . One can only observe this symmetry breaking 
at = by starting the STLS iteration with a solution 
for small but finite 9. This effectively corresponds to 
taking the limit — >■ 0, which is somewhat akin to clas- 
sical ferromagnetism, where one must consider the limit 
where magnetic field goes to zero. This stripe phase pre- 
cedes Wigner crystallization which, according to quan- 
tum Monte Carlo (QMC) calculations, occurs at [/ ~ 25 
for perpendicular fermionic dipoles in a single layer |28j . 

For 9 > arcsin(l/-\/3), the intralayer interaction devel- 
ops an attractive sliver in the plane that can eventually 
lead to collapse in the single layer P^HTT] . Here, for 
large enough U and 9, the attraction overcomes Pauli 
exclusion and the compressibility of the gas goes to zero 
ix+^i'i ~^ Oj 0) = 0)- However, we find that the bilayer 
geometry can actually stabilize the collapse to form a new 
density- wave phase that is oriented along the (j) = direc- 
tion (Fig.[T]). Referring to Fig. [2] we see that this stripe 
phase has a longer wavelength than the = 7r/2 one 
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FIG. 2: (Color online) Critical wave vector Qc/kp for the 
= 7r/2 stripe phase {0 < 9c) and the = one (6 > 6c) 
— same parameters and symbol scheme as in Fig. [l] The 
insets depict the alignment of the dipoles with the electric 
field E and the features of the two different stripe phases. 
For the = stripe phase, the density modulations in the 
two layers have a phase shift rj ~ 29, while the wave vector 
Qc decreases with increasing tilt angle 9 down to gc = for 
6 = 11/2 (filled [red] diamond), where the gas collapses. For 
density modulations along = ■7r/2, qc appears to be fixed by 
the density. 

and is dependent on geometry. Indeed, we find that qc 
smoothly decreases with increasing 9, reaching = at 
9 ~ 7r/2, where the intralayer attraction always appears 
to cause collapse at a fixed Uc- Away from 9 = 7r/2, we 
find that the <j) — Q stripe phase has qc ^ 1/d in the limit 
d — > 00, which is reminiscent of the behavior of charge 
density waves in electron- hole bilayers. 

The = stripe also features a nontrivial phase shift 
77 between the density waves in each layer. At the stripe 
transition, it can be shown that 

^^ ^ i;i2(q)[l - Gi2(q)] .„^ 
K2(q)[l-Gi2(q)]| ■ ^ ^ 

When V12 and G12 are real, like for the = 7r/2 stripe 
phase, then e*** — 1 and the density waves in each layer 
are in phase, as mentioned previously. However, V12 is 
complex for the = stripe phase and thus the density 
waves are generally shifted with respect to one another. 
Indeed, as shown below, the interlayer correlations are 
small in this phase, i.e. IG12I <C 1, therefore the phase 
shift corresponds to 77 ~ 20 (see insets of Fig. [2]) and is 
essentially independent of kpd. 

The existence of two stripe phases leads to a new quan- 
tum phase transition where the stripes change their ori- 
entation. In Fig. [1] this occurs beyond the critical point 
9c — 0.75 and Uc — 15.65 where the two stripe phase 
boundaries meet. Here, when kpd is fixed, the transition 
can be accessed by changing the tilt angle 9. Alterna- 
tively, one can fix 6* < 7r/4, which is below the onset of 
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FIG. 3: (Color online) Intra- and interlayer pair correla- 
tion functions gijiy) for increasing values of the interaction 
strength U towards the = 7r/2 stripe phase (0 = top 
panel) and the <j!> = phase (0 = 1.1 ~ O.SStt bottom panel). 



collapse in the single layer, and vary kpd^ since we ex- 
pect the critical angle to decrease with decreasing kpd. 
Eventually, at kpd ~ 1, one enters the regime where the 
physics of bosonic interlayer dimers dominates. 

Further insight into the stripe phases can be gained by 
examining the intra- and interlayer pair correlation func- 
tions gij(r) on the liquid side of the transition. For the 
= stripe phase (bottom panel of Fig. |3| , we find that 
neither pair correlation function changes significantly as 
we approach the transition. In particular, (?ii(r) only 
deviates slightly from the non-interacting case [U = 0), 
while 312(1") slowly oscillates close to one, indicating that 
interlayer correlations are small, i.e. IG12I ^ 1. This 
suggests that we can accurately model the </) = stripe 
phase using exchange correlations only. To this end, 
we construct a simplified STLS theory where we take 
G'i2(q) = and then determine the intralayer local field 
factor Gii(q) by feeding the non-interacting intralayer 
structure factor 5*0(9) ^ Io° dujll{q,u}) into Eq. 
We then evaluate the phase boundary for the (j) = stripe 
within this simplified HF theory. Referring to Figs.[l]and 
[2] we see that we obtain very good agreement with the 
full STLS calculation, particularly when U and 9 are not 
too large so that the intralayer p-wave pairing correla- 
tions are expected to be weakest [HI [T7] . In addition, 
the collapse instability at 6* = 7r/2 is unaffected by the 
other layer since the interlayer Hartree term is zero for 
q = 0. We expect one can obtain quantitatively similar 
results for the = stripe phase using the conserving 
HF approximation 34 . 



By contrast, for the (j) = tt/2 stripe phase (top panel 
of Fig. [3]), we see that correlations beyond exchange be- 
come substantial, resulting in a pronounced "correlation 
hole" for 511 (r) with increasing interaction strength, like 
in the single-layer case [TB] — note that the STLS proce- 
dure does not guarantee that gn is always positive |23| . 
and thus we sometimes obtain unphysical negative val- 
ues. The intralayer correlations also develop a substan- 
tial (f) anisotropy as we near the stripe transition. At the 
same time, the interlayer pair correlation function (712(1") 
increases at r = 0, a feature that has been ascribed to an 
imminent bound-state instability 

Indeed, the attractive part of fi2(q) always yields a 
two-body bound state composed of one fermion from each 
layer 30, .31j. Hence, any liquid phase in the phase di- 
agram contains pairing correlations and must therefore 
be superfluid (Fig. [T]) . When the size of these interlayer 
dimers Ib is smaller than the interparticle spacing, i.e. 
Ib ^ l/^F, then the system is better described in terms 
of bosonic dimers and our approach of analyzing density 
instabilities of the Fermi liquid phase is unlikely to be 
accurate. To estimate this region of phase space where 
bosonic behavior dominates, we solve the two-body prob- 
lem, E'V'k = ^■0k + / (|^wi2(k - kOV'k', where is 
the two-body wave function in terms of relative coordi- 
nates and E is the dimer binding energy. We estimate 
the dimer size as Ib ^ fi/ \/in\E\ and then determine the 
"critical" line kplB — 1 for the bosonic regime, as plotted 
m Fig. [1] (shaded region). We see that this region is well 
separated from the stripe phase boundaries and thus we 
expect our results to be reasonable for kpd — 2. How- 
ever, the presence of bosonic dimers hastens the onset 
Wigner crystallization: QMC calculations [32l l33] pre- 
dict that perpendicularly-aligned bosons will crystallize 
at f7 ~ 8. For increasing 9, the interlayer dimer be- 
comes more weakly bound until eventually the fermions 
preferentially form pairs within the same layer instead. 
With decreasing kpd, however, the regime of interlayer 
bosons expands so that it encroaches on our predicted 
stripe transitions for kpd ~ 1 and takes us beyond the 
scope of this letter. 

Our predicted stripe phases should be accessible ex- 
perimentally with cold dipolar gases. In particular, the 
bilayer distance kpd — 2 can be achieved for a typi- 
cal 2D density n ^ 1.3 x 10^ cm~^ and layer spacing 
d — 500 nm. Polar molecules such as LiCs [S] have dipo- 
lar moments D ^ 0.35 — 1.3 Debye (corresponding to 
U ^ 1 — 14) , which allows one to explore both = and 
= 7r/2 stripe phases. Furthermore, the newly explored 
NaK molecules [12] allows one to reach even larger values 
of the interaction strength (D ^ 2.7 Debye and U ~ 28). 
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